High energy emission components in the short GRB 090510 



Alessandra Corsi 1 ' 2 ' 3 , Dafne Guetta 4 and Luigi Piro 2 
alessandra. corsiOromal . infn. it 
guetta@oa-roma . inaf . it 
luigi .piroOiasf-roma. inaf . it 

ABSTRACT 

We investigate the origin of the prompt and delayed emission observed in the short 
GRB 090510. We use the broad-band data to test whether the most popular theoret- 
ical models for gamma-ray burst emission can accommodate the observations for this 
burst. We first attempt to explain the soft-to-hard spectral evolution associated with 
the delayed onset of a GeV tail with the hypothesis that the prompt burst and the 
high energy tail both originate from a single process, namely synchrotron emission from 
internal shocks. Considerations on the compactness of the source imply that the high- 
energy tail should be produced in a late-emitted shell, characterized by a Lorentz factor 
greater than the one generating the prompt burst. However, in this hypothesis, the 
predicted evolution of the synchrotron peak frequency does not agree with the observed 
soft-to-hard evolution. Given the difficulties of a single-mechanism hypothesis, we test 
two alternative double-component scenarios. In the first, the prompt burst is explained 
as synchrotron radiation from internal shocks, and the high energy emission (up to 
about 1 s following the trigger) as internal shock synchrotron-self-Compton. In the 
second scenario, in view of its long duration (~ 100 s), the high energy tail is decoupled 
from the prompt burst and has an external shock origin. In this case, we show that a 
reasonable choice of parameters does indeed exist to accommodate the optical-to-GeV 
data, provided the Lorentz factor of the shocked shell is sufficiently high. Finally, we 
attempt to explain the chromatic break observed around ~ 10 3 s with a structured jet 
model. We find that this might be a viable explanation, and that it lowers the high 
value of the burst energy derived assuming isotropy, ~ 10 53 erg, below ~ 10 49 erg, more 
compatible with the energetics from a binary merger progenitor. 
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Subject headings: gamma-ray burst: individual (GRB 090510) — X-rays: bursts - 
radiation mechanisms: non-thermal 



1. Introduction 



Traditionally divided in "long" and "s hort" on the basis of their 7-ray duration (longer or 
shorter than 2 s, iKouvelioutou et al.lll993l ). gamma-ray bursts (GRBs) are charac terized by a 



prom pt release of 7- and X-ray photo ns, followed by a multi - wavelength afte rglow (jCosta et al 



19971 ) emission. The fireball model (e.g. iMeszaros k, ReesHl992l ; ISari et al.lll998l ) explains the GRB 
electromagnetic emission as a result of shock dissipation in a relativistic flow, or "fireball" , taking 
place at distances greater than 10 -5 -10 -2 pc from the central source. Despite such a large distance, 
electromagnetic observations have revealed important clues about the progenitors, favoring two 
main models: the coalescence of a binary consisting of two neutron stars or a neutron star and a 
black hole for short GRBs; and the death of massive stars (collapsars) for long GRBs. Both of 
these are commonly believed to end in a BH-plus-torus system, where torus accretion powers the 
fireball jet. 



In the inte rnal-external shock scenario of the fireball model (see e.g. IMeszaros fc Reed 1 19931 ; 
Sari et al.1 Il998l ). GRB prompt and afterglow emissions are thought to be produced by particles 
accelerated via shocks in an ultra-relativistic outflow released during the burst explosion. While 
the prompt emission is related to shocks developing in the ejecta (internal shocks, IS), the afterglow 
arises from the forward external shock (ES) propagating into the interstellar medium (ISM). Syn- 
chrotron and synchrotron-self-Compton (SSC) emission by the accelerated electrons are typically 
invoked as the main radiation mechanisms. The Fermi satellit^] is currently bringing exciting new 
results, detecting high energy (~ GeV) extended tails whose presence is particular l y intriguing 



in th e case of short GRBs (e.g. lAbdo et al.ll2009bl : iGiuliani et al. 



2010 



Omodei 



2008; 



Ohno et al 



20091 ). Fermi observations of GRB 081024B and GRB 090510 clearly point to the existence of a 
longer-lasting high energy tail in the GeV range following the main event, motivat ing a deeper 



exploration and re-examination of the fireball physics and radiative mechanisms (e.g. lAsano et al 



2009: 


Cor si et al. 


2009; 


Gao et al. 


2009; 



Here we study the conditions under which the high energy observations of GRB 090510 can 
be accommodated within the most popular theoretical models. This work is organized as follows. 
In Sec. [2] we summarize the observations for GRB 090510. In Sec. [3] we discuss in more detail 
the spectral properties of the emission during the first second, and their implications for the com- 
pactness of the source. We test whether the complex emission observed during the first 1 s can 



1 http://fermi. gsfc.nasa.gov 
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be attributed to a single component, specifically synchrotron emission from IS. After showing that 
a single mechanism does not offer a straightforward explanation, we test an IS synchrotron plus 
SSC scenario. In Sec. 2] we examine separately the high energy emission in the context of the ES 
scenario. Finally, in Sec. [5] we give our conclusions. Hereafter we adopt as To the onset of the main 
GRB pulse which, as specified in the following section, is about ~ 0.5 s after the precursor that 
triggered the Fcrmi/GBM. Also, hereafter e e and e# are the fraction of energy going into electrons 
and magnetic fields, respectively; n is the ISM density in particles/cm 3 ; Ei so is the isotropic kinetic 
energy of the fireball; and p is the power-law index of the electron energy distribution in the shock. 



The Observations 



Characterized by a Tgo of 0.3 s (jUkwatta et al.l 120091 ). which places GRB 090510 in the short 
GRB category, the main burst was followed by an extended high energy tail, observed by both AG 



ILE/GRID and Ferm i /LAT on a timescale much longer than the prompt burst event (jAbdo et al 
2009bl : iGiuliani et all Eoiok AGILE detected GRB 090510 at T = 00 : 23 : 00.5 UT, May 10 



triggering on the sharp main peak of the GRB (about 0.5 s after the Fermi trigger on a sm aller 
precursor). The Swift/BAT also triggered around the main burst peak (jHoversten et al.ll2009l . and 
references therein) at 00 : 23 : 00.4 UT. Hereafter we adopt as To the onset of the main peak. The 
0.3-10 MeV and > 25 MeV emission of the burst showed a clear dichotomy between the low- and 
high-energy 7-ray emissions, so that two time intervals were define d: interval I from To to To + 0.2 
s, and interval II from T + 0.2 s to T + 1.2 s ([Giuliani et al.ll201oh . The main peak in the AGILE 
calorimeter (MCAL, 0.35-100 MeV) ended around 0.2 s, at which time the signal suddenly started 
to be observed in the GRID (> 100 MeV). The Swift/BAT lig ht curve showed two p ulses between 
0.2 s and 0.3 s, of amplitude much smaller than the first peak (jUkwatta et al.ll2009l ). 



The AGILE photon spectrum of Interval I is well modeled by a power law w ith index a 



— 0.651q32 and exponential cutoff E c = 2.8 MeV (see the top panel of Fig. 4 in IGiuliani et al 



2010 ). This is con sistent with the Fermi observations, which are well modeled by a Band spectrum 
([Band et al.lll993l ) with E pea k = 2. 8 MeV, a = —0-59 ± 0.04, f3 < —5, and a normalization constant 
A~8x 10" 2 ph cm" 2 s" 1 keV -1 (|Abdo et al.ll2009bh . 



During interval II, the spe ctrum undergoes a soft-to-hard evolution. Specifically, for AGILE 
(see the lower panel of Fig. 4 in IGiuliani et aliboid ). the MCAL spectrum is a power-law, N(E) ~ 



C(E /WO keVf, of photon index (3 = 
by considering that the emit t ed 0.5 MeV 



1.58i°- 13 



\l, and C = 1.6 x 10~ 2 ph cm -2 s _1 keV" 1 (derived 



erg cm 2 ; see 



10 MeV fluence during interval II was of ~ 3.1 x 10 



Giuliani et al.l (|2010i )). This spectrum is consistent with the one measured by the 
GRID, which is well fit by a power law with index f3 = — 1.4 ± 0.4, for a 25-500 MeV fluence of 
~ 2.12 x 10~ 5 erg cm" 2 (|Giuliani et al.l boioh . 



Time-resolved spectral fits of Fermi data from interval II show a progressive evolution from a 
Band plus power-law to a single power-law spectrum. Between Tq + 0.1 s and Tq + 0.3 s, the Band 
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component is still evident and has best-fit spectral indices of a = —0.48 ± 0.07, (3 = —3.09^'^. In 
the last two temporal bins the contribution of the Band component becomes less important, with 
the peak flux decreasing by about an order of magnitude. The power-law component in these bins 
has a best-fit photon index/normalization at 1 GeV of /3 = — 1.54+q'q 4 / A pow = (6. 4^^) x 10~ 9 

1 keV -1 , respectively 



ph cm s keV , and (3 



-1.92 



+0.20 



(|Abdo et alJl2009bh . 



-0.22 



pow 



(3.7; 



-1.3 > 
-1.1/ 



x 10 ph cm s 



At a redshift of z = 0.903 (jRav et all 12003 ). the 10 keV - 30 Ge V measured burst fluence 
during ~ 0.5 s since To implies an isotropic energy release of ~ 10 53 erg ( Abdo et al.ll2009bl ). which 
is extremely high for a short GRB. 

A temporal analysis of the high-energy tail at energies above 0.1 GeV as observed by the 
Fermi/LAT shows a GeV flux rising in tim e as ~ t^ bs , and decaying as ~ 



following the trigger (jGhirlanda et al.ll2010l ). Similarly, the signal detected by the AGILE/GRID 
showed that during Interval II and later on, up to 10 s after To, the emission can be described by a 
powe r-law temporal decay of index 5 = —1.30 ± 0.15 (see the top panel of Fig. 3 in lGiuliani et al 
20ld ). 



The Swift XRT began observing the burst about 100 s after the trigger (jHoversten et al.ll2009i ). 
In X-rays, a steepening is observed around 1500 s, which changes the power-law temporal decay 
index from 5 = —0.74 ± 0.03 to 5 = —2.18 ± 0.10. The o ptical light curve first ris es until ~ 1600 s, 
and then decreases as a power law with 5 = — 1.13±g;^ ()De Pasquale et al.ll2010l ). This power-law 
decay is shallower than the one observed in X-rays after about 1500 s. We also note that the optical 
emission is observed to peak much later than the extended high energy tail observed by the LAT. 



3. The first 1 s of emission 

The dichotomy of the spectral behavior observed during the first second of emission suggests 
that the properties of the source are evolving between interval I and interval II. While during interval 
II the observation of a ~ 30 GeV photon requires an optically thin source in the GeV range, during 
interval I the absence of emission above 100 MeV and the unusually steep high energy photon 
index (/3 < —5, much smaller than values typically observed in GRB prompt spectra), suggest that 
thickness due to pair production plays a role. 

The key parameter determining the optical thickness due to pair production is the Lorentz 
factor of the shell. As we show in detail in this section, the bigger the Lorentz factor, the lower 
the optical thickness. Thus, a scenario possibly explaining GRB 090510 observations could be the 
following. The central engine emits a first shell with Lorentz factor Tj, responsible for the first 
peak observed in the GRB light curve, which covers Interval I. This peak is characterized by an 
observed spectrum with no emission above 100 MeV and an extremely steep high energy spectral 
slope, thus suggesting that Tj is such that the source is optically thick above 100 MeV. Later on, 
the central engine emits a series of shells responsible for the other multiple peaks observed during 
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interval II. These shells are characterized by a Lorentz factor between Tj and Tjj, where Tjj > Tj 
and is such that the source is transparent to GeV photons. From the point of view of the physical 
properties of the source, this implies that the GRB central engine should be emitting shells with 
progressively higher velocities. The hypothesis that the source emits shells of different velocities is 
indeed the basis of the IS model. 

In the above scenario, time-resolved spectroscopy during interval II should show a progressive 
transition from an optically thick to an optically thin spectrum in the GeV range. Also, any 
spectrum obtained by integrating over multiple peaks happening during interval II would show a 
superposition of spectra emitted by shells with different T factors, progressively more transparent to 
GeV photons. Time-resolved spectroscopy by Fermi during interval II does indeed show a transition 
from a spectrum peaking around few MeV (Band component) to one with substantial emission in 
the GeV range (power-law component). Also, these components are obs erved simultaneousl y in the 



spectrum integrated from T + 0.1 s to T + 0.3 s (red curve in Fig. 2 of lAbdo et al.ll2009bl ). which 



includes at least two peaks (see Fig. 1 in lAbdo et al.l (j2009bl )) following the main one. Thus, on 



general lines, the observed spectral evolution during the first 1 s of emission from GRB 090510 is 
consistent with the hypothesis of a transition from an optically thick to an optically thin spectrum 
in the GeV range, which would naturally explain the delayed onset of the GeV tail observed by the 
LAT. 

To validate the scenario outlined above, however, a more quantitative test is necessary. It is 
worth stressing that, according to the observations, between interval I and II not only is the source 
becoming optically thin to GeV photons, but also the spectral shape of the observed emission 
is changing substantially. In particular, a crucial point to verify is whether the evolution of the 
Lorenz factor required to justify a transition toward a smaller thickness in the GeV range also 
agrees, within the IS model, with a shift in the spectral peak from a few MeV to more than 1 GeV 
(as observed in the transition from a Band to a power-law spectrum). In what follows, we analyze 
this scenario in detail. 



3.1. Interval I 



3.1.1. Thickness to pair production 



We can argue that during interval I, and in particular between Tq and Tq+0.1 s (see lAbdo et al 



2009b! ) . the unusually steep high-energy photon index observed by Fermi (/3 < —5) is due to optical 
thickness from pair production on an underlying Band function with A, E pea k and a as the observed 
ones (see Sec. [2]), but with /3 > — 5. More specifically, hereafter we make the hypothesis that the 
true (unabsorbed) high energy spectral index has a value of = —3. 5. In fact, according to the 
complete spectral catalog of BATSE bright GRBs (jKaneko et al.ll2006l ). the tail of the distribution 
for the spectral index of well-modeled spectra is around f3 = —3.5. While the BATSE catalog did 
not show any significant difference between the spectral parameters of short GRBs and those of 
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long ones, we note that two of the short GRBs in that sample indeed had (3 ~ —3.5 (see Table 14 in 
Kaneko et al.ll2006l ). Moreover, being in the tail of the distribution, f3 = —3.5 would reconcile GRB 
090510 observations with the more commonly observed properties of GRB prompt spectra, while 
minimizing the implied value of the r 77 for pair production. In fact, we cannot have r 77 >> 1 if 
the observed spectrum is non-thermal and the light curve shows high temporal variability. We also 
note that an unabsorbed j3 of —3.5 woul d be consistent, wi thin the errors, with j3 = —3.09 



+0.21 
-0.35 



Of 

the Band component observed by Fermi (jAbdo et al.ll2009bl ) during To + 0.1 s and To + 0.2 when, as 
discussed in Sec.El we expect a contribution from less thick shells whose spectrum (in the observed 
energy band) is evolving from a Band to a power-law shape. This of course should be taken with 
the caveat that even j3 = — 3.09^q 35 could still be affected by absorption. 



The high-energy part of a Band spectrum reads 



N(E) = 0.08 



(a - (3)E peak 
e(2 + a)100 keV 



E 



ph 



100 keV ) cm 2 s keV 



CBand 



E 



100 keV J cm 2 s keV 



ph 



(1) 

with a ~ -0.59, E ppnk ~ 2.8 MeV for GRB 090510 (see Sec. EJ). T his equation is obtained from 



Eq. (1) of Band et all (|1993l ) by using E penk = (2 + a) En (see e.g. |Piranlll999h . Note also that 
the multiplicative factor e@~ a in Eq. (1) of iBand et al.l ( 19931 ) is included in the first factor in 
parenthesis in our above equation. If the true spectrum has f3 = —3.5, then 



C 



Band 



-0.59+3.5 



and we should have 



-0.59 + 3.5)2.8 MeVY 
e(2 - 0.59)100 keV / 



r 77 (100 MeV) > 4 



(2) 



(3) 



to reco ncile this spectrum wit h the observed j3 < — 5. The r 77 for pair production is expressed as 
follows jLithwick fc Sarill200lh : 



0.1a T N 7>Ean{E) 
4ttR 2 



(4) 



In the above relation, <jt is the Thompson cross section, R is the size of the source, and N^ >Ean ^ 
is the number of target photons; i.e., the number of photons with energy above E an , where 



Ean(E) 



(Fm e c 



2\2 



2.6 x io 5 r 2 



keV. 



(5) 



E(l + z) 2 (£/keV)(l + z) 2 

This accounts for the fact that a photon with energy E in the observer frame may be attenuated by 
pair production through interaction with softer photons, whose energy (also in the observer frame) 
is equal to or greater than E an {E). Also, for a power- law spectrum of the form 

ph 



one has 



N. 



~(>E an (E) 



N{E) = C(E/100 keV f- ,, , 

cnrs keV 

C4vr((i L /cm) 2 (^/s)(TJ QW (E)/keV) 1 +/ 3 
-(l + ^XlOO^l + z) 2 



(6) 



(7) 



-7- 



where we are supposing (5 < — 1. For convenience, we define E max as 

T-yyiEmax) = 1 . (8) 

It is evident from Eqs. and (|7|) that r 77 scales with energy as t a E~^ 1Jr ^; the requirement 
r 77 (100 MeV) > 4 (see Eq. (|3|)) then allows us to constrain the value of E max as follows: 

1/4 > (E max /W0 MeV)-^- 3 - 5 ) =► E max < 4' 1 / 2 - 5 x 100 MeV ~ 60 MeV . (9) 

This requirement on E max implies the following condition on the Lorentz factor of the shell. Using 
R = 2cT 2 5t obs /(l + z) = 6 x 10 10 r 2 [5t obs /((l + z)s)] cm and substituting Eqs. © and © into 
Eq. we have 



1.8 x 10- 47 C(d L /cm) 2 (2.6 x 10 5 ) 



(E max )(l + z)(2+2«(100WWs)(-l " P){E max /keV)W) 



T. 



1/(2-2/3) 



(10) 



28 



For interval I, setting C = Csand = 584, Tjj(E max ) = 1, (3 = —3.5, z = 0.903 and = 1.8 x 10 
cm as appropriate for GRB 090510, we get 

E max < 60 MeV, (11) 
T 7 ~ 160(<5i o fo s /100 ms)~ 1 / 9 (£' max /60 MeV) 5 / 18 . (12) 



3.1.2. Thickness for scattering on pairs 



To explain a non-thermal spectrum and a high temporal variability, the number of pairs created 
by the optically thick portion of the spectrum should remain small. This is in order to avoid the 
Thompson op tical depth for photon scattering on the created pairs becom i ng much greater than 
unity (see e.g . 



Sari &: Piran 



2009a 



Pe'er &: Waxman 



2004 



Guetta et al. 



2001 



Lithwick fc Sari 



2001 



Abdo et al 

19971 ) . Since we can reasonably assume that each photon of E > E max creates a pair, 



the number of pairs is approximately 



Npair ~ 



N. 



~/>E„ 



The Thompson optical depth is thus (|Abdo et al.ll2009ar ) 



T^±{E m ax) ~ 



4vri? 2 



(13) 



(14) 



We note that the above expression for r 7 ± was also used bv lAbdo et al.l (|2009al ). who pointed 
out that it's typically difficult to have a source optically thin to scattering on pairs when the optical 
thickness to pair production is very high. Using Eqs. ([7]) and ()14|) . we can write 

a T C (dL/cmfiE^/keY) 1 ^ 



r 7 ±(^„ 



4 C 2(-l-/3)(100)' 3 ( ( 5Ws)r4 



2 x 10 9 (<5W100 ms)- 1 ( J B max /60 MeV)" 2 -^" 



where we have used (5 
thus implies 



(15) 

-3.5, z = 0.903, C = C Ba nd = 584, d L = 1.8 x 10 28 . Requiring r 7± < 1 



1/ > T 7± = 200(« ofa /100 ms)- 1 / 4 (S max /60 MeV)" 



•5/8 



(16) 
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3.1.3. Thickness for scattering on electrons 



Further constraints on the Lorentz factor come from scattering of the emitted photons on 
electrons inside the shell. According to the IS model, a fraction e e of the internal energy of shocked 
particles goes into accelerating the electrons. The shock- accelerated electrons then radiate via 
synchrotron (and IC) emission. A necessary condition for radiation from IS to be observed is 
that the source is optically thin for photon scattering on electrons associated with baryons present 
inside the shell itself. When the r 7e for photon scattering on electrons is high, the spectrum of the 
observed radiation is modified by the standard assumptions of thin synchrotron and IC emission, 
and effects related to the presence of the so-call ed electron photospher e need to be considered . 
Another condition we thus need to set is that (IRees k. Meszarosl 120051 ; IPe'er Waxmanl |200J; 
Guetta et aDhoOll : iLithwick fc SarilhoOll : iMeszaros & Reeslboooh 

on 



where we have indicated with 



N, 



AttR 2 

L5t b s 



<1, 



baryon 



(17) 



(18) 



(1 + z)e e Tm p c 2 

the number of electrons associated with baryons inside the shell, and R ~ 2cT 2 5t b s / (1 + z) is the 
radius of the shell. We thus get 



r 7e ~ 10 9 L 52 (1 + zXJWlOOms)-^-^- 1 < 1 • 



(19) 



By integrating in the 10 keV - 100 MeV energy range a Band function wit h normalization 
constant ~ 0.08 ph/cm 2 /s/keV, E peak = 2.8 MeV, a = -0.59 (jAbdo et all l2009bh . /3 = -3.5, and 
multiplying by Aird 2 ^ (with z = 1.903), we estimate the luminosity in the cosmological rest frame 
to be L ~ 3 x 10 53 erg/s. Thus, 



r/>r 



7e 



140 e- 1/5 (St obs /100 ms) 



-1/5 



(20) 



3.1.4- Synchrotron emission from IS 



As we have seen in the previous sections, the spectrum observed during interval I could be 
reconciled with more commonly observed Band spectra by requiring the optical thickness for pair 
production to be responsible for the unusually steep high-energy spectral decay. Here we analyze 
the conditions under which such a spectrum could be expla ined as synchrotron em ission from IS. 
The peak of the synchrotron component in the IS model is (jGuetta Sz Granotll2003l ) 

t:,„ „/. ~ I :l [ ) el^e^L^T^idto^/im ms)" 1 GeV . (21) 



p-1 



For the case of GRB 090510, a high energy slope of /3 = —3.5 can be explained in the (optically 
thin) IS scenario by requiring p ~ 5 (so that the expected photon spectral index is — p/2 — 1 = —3.5; 
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see iGuetta fc Granotl 12003 ) and by requiring that Eq. (|12p is valid. Using these, and setting 



€b ~ 0.5, L52 = 30, we can write Eq. (|2ip as follows: 

E P eak = 0.3(5t obs /W0 ms)- 7 l\E max /W MeV)- 5 / 9 MeV . 



(22) 



To have a peak around 1 MeV, we need St b s ~ 20 ms and E max ~ 60 MeV. For such a 
value of the variability timescale, from Eq. (fT2l) we get Tj ~ 200, and using Eqs. ([TBI) and ([20]) 
we also have r 7 ± ~ 300 and T 7e ~ 200. This implies that r 77 (60 MeV) ~ 1 ~ r 7e , but r 7 ± ~ 6. 
We thus expect to have some effects from scattering of the emitted photons on the created pairs. 
Numerical simulations are the best way to predict these effects, since there are different processes 
that come into play during the dynamical timescale. In addition to pair production and scattering 
on electrons, discussed before, one s hould also consider, e .g., re-heating of the electron population 
due to synchrotron self- absorption (|Ghisellini et al.lll988l ). and pair annihilation. All these effects 
combined together can modify the observed spectrum. 



Pe'er &: Waxmanl (J2004J) have carried out time-dependent numerical simulations within the IS 



model, describing cyclo-synchrotron emission and absorption, inverse and direct Compton scatter- 
ing, and pair production and annihilation (including the evolution of high energy electromagnetic 
cascades), allowing a calculation of the spectra resulting when the scattering optical depth due to 
pairs is high, thus prese nting deviations from the simple p r edicti ons of the thin case IS model (e.g. 
Guetta &: Granotl l2003l ). In particular, IPe'er &: Waxmanl (|2004l ) have shown that from moderate 
to large values of t 77 , the resulting spectrum peaks in the MeV range (as was the case for GRB 
090510), shows steep slopes at lower energies, and exhibits a sharp cutoff at ~ 10 MeV. For large 
compactness, scattering by pairs becomes the dominant emission mechanism, as we have seen here 
for GRB 090510 (r 7 ± ~ 6). In such a case, electrons and positrons lose their e nergy much faster 
than the dynamical timescale, and a quasi-Maxwellian distribution is formed (IPe'er &; Waxman 
2004). The energy gain of the low-energy electrons by direct Compton scattering results in a spec- 
trum steeper than Maxwellian at the low-ener gy end, indicating that a steady state did not develop. 
Predicted slopes in vF v are 0.5<2 + a<l ([Pe'er fc Waxmanl 120041 ) . 



In the case of GRB 090510, the low energy photon spectral slope is a = —0.59. This value 
implies a very steep rise in the vF v spectrum, of about ~ 1.4. Detailed modeling of the spectrum 
for high compactness is beyond the purpose of this paper. These considerations, however, allow us 
to conclude that, overall, the spectrum observed during Interval I may be accommodated within a 
high compactness, synchrotron IS scenario. 



3.2. Interval II 



3.2.1. Transparency to GeV photons 



During interval II, photons up to 30 GeV were observed by the Fermi/LAT (jAbdo et al. I l2009bh . 
Thus, in contrast to interval I, the source should be optically thin to pair production and have 
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r 77 < 1 at E max > 10 GeV. During Interval II, the observed spectrum is consistent with (see 
Sec. ED 

N(E) = C pow (E/W0 keVf ph - , (23) 

cm z s kev 

where /3 1.58 and C pow ~ 1.6 x 10" 2 . Thus, using these values in Eq. (fTUj) . and setting z = 0.903, 

d L = 1.8 x 10 28 CHI, 7~/y/y (E max ) = 1, and E max = 10 GeV, we obtain the following requirements: 

E max > 10 GeV; (24) 
Fjj ~ 430(5^/100 ms)- 1 / 5 - 16 ( J E maa; /10GeV) - 58 / 5 - 16 . (25) 



3.2.2. Synchrotron emission from IS? 

If we assume that the spectrum observed during Interval II is dominated by synchrotron 
emission from IS; i.e., it is generated by the same radiation mechanism explaining the emission 
in Interval I, then it is necessary to require that the peak of the synchrotron component has a 
soft-to-hard evolution, with E peak ~ 2.8 MeV during interval I and E peak > 1 GeV during interval 
II. Substituting Eq. ([25]) into fl2T} we have 



E peak ~ 3 x 10- 5 (5W100 ms)-°- 61 (£ ma:E /10 GeV)" ' 22 GeV , (26) 

where we have set e e ~ es ~ 0.5, p ~ 5. We have also used the fact that during interval II, most 
of the emitted energy is in the GRID ener gy range, with a measured 25 MeV - 500 MeV fluence of 



2.12 x 10~ 5 erg cm" 2 ([Giuliani et al.ll2010h . thus giving L ~ 4^(2.12 x 10~ 5 erg cm" 2 / 1 s ) ~ 1q53 



erg/s during interval II. From the above equation, it is evident that even setting St Q b s ~ 1 ms, we 
have E peak « 1 GeV for E max > 10 GeV. 



3.2.3. SSC emission from IS: a better explanation 

The extreme soft-to-hard evolution observed between interval I and II suggests an alternative 
two-component explanation. Specifically, one could think of the emission in interval I being domi- 
nated by IS synchrotron of a slower shell with r ~ Tj (see Sec. 13. If) . while the emission in interval 
II dominated by IS SSC of a late-emitted faster shell (r ~ Tjj > T/), whose SSC component falls 
in the observed band, while the synchrotron counterpart is shifted to lower energies. 

In Fig. Q] we show a possible solution within this scenario: during interval II, the high energy 
emission is dominated by the IC component of a faster shell with Tjj ~ 645, L§2,ii ~ 6, 5t b s ji ~ 15 
ms, e e ji = 0.35, €b,ii = 0.008, pjj = 4.8. We stress that what we show in this figure implies 
that a viable parameter choice does exist to accommodate the observations within this model. 
However, such a solution is not necessarily unique, and a larger parameter range may exist. For 
a value of p = 4.8, we expect a high-energy spectral slope of /3 = — 1 — p/2 ~ —3.4 for the 
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synchrotron component photon spectrum, consistent with our initial hypothesis that the true high- 
energy spectral slope is /3 = —3.5, and it's initially (between To and To + 0.1 s) made steeper 
(/3 < —5) by opacity due to pair production . We note that the slope observed by Fermi between 
T + 0.1 s and T + 0.3 s (/3 = -3.09+" f 5 , kbdo et al.ll2009bh would agree with the hypothesis 



that we expect the spectrum to become progressively more transparent (see Sec. 13.1.1]) . We should, 
however, keep in mind that between To + 0.1 s and To + 0.3 s some effects due to absorption might 
still be present. According to the predictions of the IS model, we also expect a photon index of 
— 1.5 for the SSC component, which a grees with the value of /3 = -1.581^? observed by AGILE 
during interval II (jGiuliani et al.ll2010i ). 



We finally underline that other interesting scenarios have been propose d to explain the high- 
energy emission observed during interval II. For example, iToma et al.1 (120101 ) recently showed that, 
in the framework of the IS model, effects related to up-scattered photospheric photons may become 
visible, and explain the delayed high-energy tails observed by the Fermi/LAT. The explanation we 
propose here is thus limited to considerations based on the (simpler) assumption of an optically 
thin IS model. But, indeed, other scenarios are possible. 



4. Synchrotron emission from the ES 



In this section, we test whether the high-energy emission observed by the Fermi/LAT, and the 
optical-to- X-ray emission observed later on by Swift, can be explained as ES afterglow, while the 
emission in the Fermi/GBM, AGILE/MCAL and Swift/BAT is due to IS. In this way, one can easily 
account for both the high temporal variability observed during the prompt burst (as related to IS), 
and for the delayed onset of the high-energy emis s ion (a s related to the onset of the afterglow). 
This hypothesis, first proposed by lGhirlanda et al.1 pOlCj) on the basis of the temporal behavior of 
the high energ y tail observed in the LAT up to 100 s after the burst, was then confirmed as a viable 
possibility by iDe Pasquale et al.l (120101 ) performing a broad-band analysis based on Swift BAT, 
XRT, UVOT, Fermi GBM, and LAT data. The striking feature of the broad-band observations is 
that the spectral ene rgy distribution of the e mission observed at 100 s is consistent with a single 
spectral component ( De Pasquale et al. 201ol ). Here we assume the most natural hypothesis of it 
being simply the synchrotron high energy tail. 



De Pasquale et al.l (120101 ) have suggested that, within the ES model, the peak observed around 



~ 0.2 — 0.3 s after the BAT trigger in the Fermi/LAT light curve of the extended tail could be 
associated with the fireball deceleration time, while the peak observed in the optical range could be 
due to the synchrotron peak frequency v m crossing the band. In light of these considerations, we 
have modeled the synchrotron emission from the ES to test if a reasonable set of parameters does 



Alternatively, one could suppose that all the optical-to-GeV emission is generated by SSC of a synchrotron IS or 
ES component peaking at much lower energies, but this would be a rather non-standard scenario, which we do not 
analyze here. 



-12- 



indeed exist to provide such an explanation. To this end, we adopt the prescriptions by I Sari et al 
(|1998l ) for the peak flux f m , the injection frequency v m , and the cooling frequency u c : 



-2 . 



(27) 



Vm. OC 



p — 1 



-i 



(28) 



(29) 



We further rescale the expression for v c by a factor of Y 2 , to account for the effect of SSC 
l osses on the syn chrotron spectrum. Here Y is the Compton parameter and is defined as follows 
rtSari fc EsinlboOlh : 



Y 



Y 



max 



max 




1, 



1,1^ 



in fast cooling ; 



l/(4-p) 



-(p-2)/2(4-p)- 



in slow cooling , 



(30) 
(31) 



where v c is the non-rescaled value of the cooling frequency, as in Sari et al. (jl998l ) and Eq. ([29]) . To 
model the behavior of the high-energy tail, we consider the whole evolution of the Lorentz factor 
r of the shell, using an approximate sharp transition from the coasting phase, when 



r 



to the deceleration phase, when (e.g. I Sari et al .11 19981 ) 

r(£ fts) = r (t o bs/tdec) 3 ^ 8 



Here td ec is the deceleration time in the observer's frame, given by (pari fc Piranlll999l ) 



tde 



32vrrgnm 



1/3 



(1 + ^), 



(32) 



(33) 



(34) 



where m„ is the proton mass. 



In Fig. [2] we show what we obtain for the parameter choice Tq 
3 x IP" 3 , e P = 0.1, E ip n 



10 4 



n 



10" 



3.7 x 10 53 erg. These values are consistent with the results by 
Kumar fc Barniol Duranl (|2009l ). in which the GeV light curve was modeled starting at ~ 1 s 



after the BAT trigger, while here we are modeling also its peak around 0.2-0.3 s. The high value 
of the Lorentz factor is required to have i^ec ~ 0.3 s, so as to explain the peak observed in the 
LAT light curve. We also note that while the very low density value is still consistent with those 
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that can be expected around short GRBs in the coalescing binary progenitor scenario (see e.g. 
Belczvnski et al.ll2006l ). the isotropic energ y is much higher (th ough comparable to the one derived 
from the fluence observed in the LAT, see lAbdo et al.l (|2009bl )). 



In the hypothesis that the steepening observed in the XRT light curve is due to a jet break, we 
have further evolved the Lorentz factor as T oc (tobs/tj)' 1 ^ 2 (see e.g. IPeng et alJhoOsl. and Fig. 3), 
finding that tj ~ 1000 s in good agreement with the data. Using the relation (e.g. lRhoadslll997l ) 



rfo) = 1/0; , 

and considering Eq. (|33l) . we can constrain the jet opening angle to be 

Oj = To 1 (lti i s/0.3 s) 3/8 ~0.1°. 



The energy in the jet is thus 



Ej = (e 2 J2)E iso ~ 7 x 10 47 erg, 



(35) 



(36) 



(37) 



which is more easily explained in a binary merger model. 



We note, howe ver, that after the X-ray b reak the optical flux decreases with a slope shallower 
than the X-ray one. iDe Pasquale et al.l (|2010l ) and lKumar <fe Barniol Duranl (|2009l ) have suggested 
that this may be explained by a jet break made shallower from the passage of v m through the 
optical band. With our choice of parameters, v m is crossing the optical band around the jet break 
time, and the light curve decay is still too steep (at least using our simple approximation of the 
Lorentz factor evolution). We therefore test the alternative hypothesis of a two-component jet, with 
a narrow jet component explaining the early time emission, and a wider comp onent contr i buting 
at late times to explain the excess observed in the optical band. For example, IPeng et al.l ((2005) 
considered such a model to explain t he optical light cur ve of GRB 030329. By assuming for both 
jet components the same e e , eb and p. IPeng et al.1 (|2005l ) found that the addition of a wider, slower 
component with Tq w ~ (l/10)ro n and Ej w ~ 4i?, n could explain the late-time optical excess 



observed in the lig ht curve. A structured jet model has also been invoked in other cases (e.g. 
Racusin et al.ll2008l ) to explain chromatic jet breaks. 



In light of these considerations, we have attempted to explain the optical excess observed in 
the case of GRB 090510 after t b s ~ 10 3 s by adding the contribution of a wider jet component. We 
find that the choice E Wjiso = 1.5 x 10 53 , F , w ~ (l/45)r ,„ = 220, 6 jjW = 0.4°, and p = 2.5 (with 
the other parameters left unchanged) can account for the excess observed in the optical, with little 
contribution in the X-rays (see Fig. 2 and Fig. 3). This choice also implies Ej >w ~ §Ej n . The 
chosen value of p = 2.5 is larger than the one adopted for the narrow component. This is motivated 
by the fact that the X-ray decay observed before 10 3 s, dominated by the narrow component, is 
shallower than the one observed after 10 3 s in the optical band (5 = —1.13^'^), which we model 
as the emission from the wider component. With p = 2.5, for v m < v op t < v c one gets a predicted 
value of the temporal decay index of — 3/4(p — 1) ~ —1.12, in agreement with the observed one 
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within the uncertainties. Incidentally, we note that for the case of GRB 080319B, iRacusin et al 
(|2008l ) obtained different p values for the narrow- and wide-jet components, as we are finding here 



We finally test whether, for our choice of parameters, the contribution of SSC emission to 
the observed flux is indeed ne gligible (as suggested by the SED at 100 s being consistent with a 
single spectral component; see lDe Pasquale et al.1 (J201CJ)). The peak flux of the SSC co mponent, in 
the T homson limit, is related to the synchrotron one by /^p ~ 10~ 6 n(i?/10 19 cm)/ m ( Sari Esin 



20011 ). In our case, the peak of the synchrotron component is constrained to fit the optical flux 
measured by the UVOT, which is about 10 Jy (see Fig. [2]). This means that for n = 10 -6 , the 
SSC component has a flux level below 10~ 16 Jy at all energies throughout the evolution, so that 
its contribution to the light curves plotted in Fig. [2] is completely negligible. 



5. Conclusion 

We have analyzed GRB 090510 in the context of the synchrotron IS and ES scenarios. We 
first attempted to explain the soft-to-hard spectral evolution associated to the delayed onset of a 
GeV tail with the hypothesis that both the prompt burst and the high-energy tail originate from 
synchrotron emission of electrons accelerated by IS. Considerations of the compactness of the source 
lead us to conclude that the high-energy tail should be produced in IS developing in a late-emitted 
shell, characterized by a Lorentz factor of the order of T ~ 700, greater than the one generating 
the prompt burst (r ~ 200). However, this condition on the Lorentz factor implies a hard-to-soft 
evolution of the peak frequency of the IS synchrotron component, which does not agree with the 
observed soft-to-hard evolution. 

Given the difficulties of explaining the prompt and delayed high energy emission with a single 
mechanism (synchrotron emission from IS), we then tested two double-component scenarios. In the 
first, the emission observed during interval I is explained as synchrotron emission from IS, while the 
high-energy tail observed in interval II is explained as SSC emission from IS. In the second scenario, 
the high energy emission observed by the LAT is decoupled from the prompt burst, and has an 
external shock origin. This last scenario has the advantage of explaining in a simple way the smooth 
temporal behavior of the high-energy tail, up to 100 s after the burst, and the consistency of the 
broad-band SED observed at 100 s with a single spectral component. In the ES scenario, we show 
that a reasonable set of parameters does indeed exist to explain the optical-to-GeV observations of 
this burst, despite a high Lorentz factor being required to have the fireball entering the deceleration 
phase as early as t b s ~ 0.3 s (when the emission in the LAT is observed to peak). The ES scenario 
thus seems to account more naturally for the observations, even if a more detailed modeling of the 
late-time chromatic break is required. We have suggested that a structured jet may indeed be a 
viable explanation of such a chromatic feature. 

In conclusion, we stress that the high Lorentz factor implied by the ES scenario has some 
relevant consequences in relation to the physics of the central engine. The commonly accepted 
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fireball model invokes a series of shells expanding outward with Lorentz factor of the order of a few 
hundred, where IS first generate the prompt emission, and then, with the merged shell continuing 
to expand outward toward the external medium, an ES generates the afterglow. If, on the other 
hand, the high-energy tail is attributed to the ES emission, this implies that the source emits first 
a very fast shell, which impacts on the external medium creating an early afterglow, plus a series 
of slower shells that catch up with each other generating the prompt 7-ray emission. At the end 
of the IS phase, the merged, slower shell (with a more typical Lorentz factor of the order of a 
few hundred), would also decelerate and eventually generate an afterglow by interaction with the 
external medium. In this respect, we note that the additional component required to explain the 
shallow decay observed at late times in the optical band could be related to the ES generated by 
such a slower shell. 

The authors are grateful to Patrick Sutton for helping improve the manuscript by carefully 
proofreading its final version. A. C. thanks the Italian L'Oreal-UNESCO program "For Women 
in Science", and EGO- European Gravitational Wave Observatory, for support. The authors also 
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L„=6 e e =0.35 e B =0.008 p=4.8 r = 645 t vobs =15ms 
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Fig. 1. — Mode l ing o f the high energy emission during interval II in the IS scenario 
(buetta fc Granot) (2003) . The red arrow marks the level of BAT upper-limits (see e.g. 



De Pasquale et al.l I201Q ) while t he green diamond marks the level of the flux observed by the 



LAT around 1 GeV (see Fig. 2 in lAbdo et alJl2009bl ). 
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Fig. 2. — GRB 090510 bro ad-band modeling i n the synchrotron ES scenario. Data are taken from 
the light curves in Fig. 1 of lDe Pasquale et al.1 (|2010i ) , where the mean flux measured in each energy 
range for the different instruments (Fermi/LAT, Swift/XRT and Swift/UVOT) has been rescaled 
to give the measured flux at each specific frequency (1 GeV, 1 keV and 10 15 H z). This is done by 



requir ing at 100 s the specific flux value reported in the SED plotted in Fig. 2 of lDe Pasquale et al 



(|2010l ). The black, red, and purple solid lines represent the model predictions at 1 GeV, 1 keV and 
10 15 Hz, respectively. Parameters are set as follows, for the narrow and wide jet components (dash- 
dotted lines), respectively: e e = 0.1, e B = 3 x 10" 3 , n = 10" 6 , r , n = 10 4 , E i80 ,n = 3.7 x 10 53 erg, 
p n = 2.3, 0j, n = 0.12°, T , w = 220, E is0jW = 1.5 x 10 53 erg, p w = 2.5, 9 jjW = 0.43°. The blue lines 
represent the contribution of the narrow and wide jet components at the middle of the BAT energy 



band: throughout the evolution this is below the data/upper-limits reported in lDe Pasquale et al 



(|2010l ). in agreement with our hypothesis that the emission observed in the GBM/BAT/MCAL 
energy range should be due to IS rather than to ES. 
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Fig. 3. — Temporal evolution (black lines) of the Lorentz factor for the narrow (left panel) and 
wide (right panel) jet components. The red-dotted lines are plotted for comparison, and correspond 
to the standard evolution T oc i~ 3//8 for an adiabatic fireball expanding in a uniform medium (e.g. 



Sari et al.l Il998l ). The initial values Tq are set to be 10 4 and 220, for the narrow and wide jet 



components, respectively. After tj, which is set to be ~ 1000 s for the narrow jet (corresponding 
to 9j ~ 0.1°), and ~ 2 x 10 4 s for the wide jet (corres ponding to 6j ~ 0.4°), the Lorentz factor is 
evolved following a temporal scaling of T oc t~ 1 ^ 2 (e.g. 



Peng et all 12003 ) . 



